Theory of Diffusion Controlled Growth 
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We expand upon a new theoretical framework for Diffusion Limited Aggregation and associated 
Dielectric Breakdown Models in two dimensions [R. C. Ball and E. Somfai, Phys. Rev. Lett. 89, 
135503 (2002)]. Key steps are understanding how these models interrelate when the ultra-violet 
cut-off strategy is changed, the analogy with turbulence and the use of logarithmic field variables. 
Within the simplest, Gaussian, truncation of mode-mode coupling, all properties can be calculated. 
The agreement with prior knowledge from simulations is encouraging, and a new superuniversality 
of the tip scaling exponent is discussed. We find angular resonances relatable to the cone angle 
theory, and we are led to predict a new Screening Transition in the DBM at large 77. 
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I. INTRODUCTION 

Diffusion controlled growth first attracted attention in 
the literature on solidification, where the advance of a 
solidification front can be limited by diffusion of either 
latent heat or compositional excess ahead of the front. 
Under these conditions a planar front is linearly unstable 
with respect to long wavelength corrugation, the MuUins- 
Sekerka instability leading to a rich variety of prob- 
lems in pattern formation. Viscous fingering, arising 
when a viscous fluid is driven through a porous medium 
by less viscous one, is recognised as being a problem in 
the same class. The Diffusion Limited Aggregation model 
(DLA) of a rigid cluster growing by the irreversible ac- 
cretion of dilute diffusing particles, introduced by Witten 
and Sander Q focussed attention on the extreme limit 
of these problems, where all of the diffusion is ahead of 
the growth and quasi-static, with the added simplifica- 
tion that the MuUins-Sekerka instability applies on all 
lengthscales above the size of the accreting particles. 

Mathematically these problems share the same general 
form for the equations governing their growth, with their 
local interfacial velocity controlled by a conserved gradi- 
ent flux, 



Vn OC |9„(/)|'' , 



0, 



^interface 



(1) 



where [at least naively] 77 = 1 jsj. The generalisation 
to a range of positive 77 was introduced by Niemeyer, 
Pietronero and Wiesmann Q to model dielectric break- 
down patterns, and in our recent letter we introduced 
the idea that this can support equivalences between mod- 
els where the MuUins-Sekerka instability is controlled lo- 
cally in mathematically quite different ways. 

The DLA model has attracted enormous attention be- 
cause it contains no limiting lengthscale (save the par- 
ticle size) and so pattern formation must continue non- 
trivially on all larger lengthscales, the MuUins-Sekerka in- 
stability ruling out simple planar growth. Theoretical in- 
terest has been fueled by the fractal and multifractal |^ 
scaling properties of the clusters produced, with contro- 
versial claims [|[ ^ (and counter-claims ]|ll|, 13 ) 



an overall theoretical framework to understand the prob- 
lem. A simple mean field theory docs not capture 
the fractal aspects, which are better understood through 
various relations between exponents [|l^, ^ [T^ and [Q . 
The Cone Angle Theory gives a plausible argument 
for the fractal dimension of DLA, whilst the Screened 
Growth Model |1^ and Makarov's Theorem |Q give in- 
sight into how the multifractal spectrum of the growth is 
generated. 

The presence of a cut-off lengthscale a below which the 
physics dictates smooth growth is a crucial ingredient of 
diffusion controlled growth; it is known that otherwise 
infinitely sharp cusps develop in the interface within fi- 
nite time [^. In DLA this cutoff is fixed and set by 
the size of accreting particles, but in solidification and 
also in viscous fingering it is set by surface free energy 
leading to a local offset in the value of ^interface: this is 
dependent on the local interfacial curvature and leads to 
a velocity-dependent selection of approximately constant 
va^ for growing tips ||2^. We will gather these differ- 
ent possibilities together with the generalised cut-off law, 
that advancing tips have radius 



(2) 



of anomalous scaling, and by the long-standing absence of 



In terms of m, simple DLA corresponds to to = and 
solidification and viscous fingering correspond to m = 
1/2; in the theory below in two dimensions we will map 
onto the case where a is such that each growing tip has 
fixed integrated flux, corresponding to m = 1. 

This paper explains and expands the new theory an- 
nounced in our recent letter ||^. In Sec. || we estab- 
lish mappings between models with different 77 and to. 
This opens up the opportunity to discuss the full class 
of models via the section at m = 1, which we will show 
in Sec. |^ is particularly amenable to continuum theo- 
retical description in two dimensions. Section III dis- 
cusses the other key input to Sec. |^ that we regard 
diffusion controlled growth as a turbulence problem with 
self-organising fluctuations. In Sec. ^ we discuss how 
to cast the theory of Sec. IV in renormalised form, with 
divergent factors factored out consistently. Closure ap- 
proximations are required to obtain explicit theoretical 
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77, w growth rule 77', m' growth rule 

FIG. 1; For equivalence between growth models with different 
values of rj and m, we consider a pair of realisations matched 
down to but not including tip radii. The equivalence then re- 
quires that the two models agree about the relative velocities 
of competing tips. 



predictions, and in Sec. VI we show how the simplest 



Gaussian closure leads to a complete theory of the frac- 
tal and multifractal scaling. This turns out to be quan- 
titatively quite accurate for the zone of active growth. 
In Sec. VI] we show how one outstanding exponent, the 



tip scaling of the harmonic measure, can be pinned down 
through the use of the Electrostatic Scaling Law leading 
to a very surprising predict ion w ith which numerical data 
seems compatible. Section VIII gives a more detailed dis- 
cussion of what happens to the DBM exponents at large 
77. In Sec. IX we show examples of how the theory can be 
deployed to tackle deeper quantities such as the relative 
penetration depth in DLA, which has been the subject of 
several numerical studies and some controversy. Angular 
resonances appear leading us to draw parallels in Sec. XI 
with the earlier Cone Angle Theory. 



II. SCALING PROPERTIES AND MAPPINGS 
BETWEEN MODELS AT DIFFERENT 77 AND m 

It is central to fractal (and multifractal) behaviour in 
DLA that the measure given by the diffusion flux [den- 
sity] j = dn(j} onto the interface has singularities |^], such 
that the integrated flux onto the growth within distance 
r of a singular point is given by 



~ {r/RY 



(3) 



where R is the overall linear size of the growth. Multifrac- 
tal scaling of the flirx density leads to a whole spectrum 
of a values, with the number of regions of lengthscale 
r having /^(r) ^ {r/R)°' varying as {r / R)^^"\ but in 
the following we focus particularly on advancing tips and 
their associated exponent value atip. 

Applying this phenomenology to the scaling around 
growing tips, we can establish an equivalence between 
models at different rj and m by requiring that the rela- 
tive advance rates of different growing tips are matched. 
Consider two growths, growing governed by different pa- 
rameters {rj, m) and (77', m') respectively, which at a given 




FIG. 2: According to the mapping (^) models lying along 
any given line shown are fundamentally related and should 
have equivalent scaling properties. They should therefore be 
classifiable in terms of the value rim corresponding to any 
chosen reference value of m, such as rio or 771. For simplicity 
the graphic has been plotted for d — 2 and taking a = 2/3 
completely independent of 77 (see later): if a does vary from 
line to line, then the lines will not be confocal. 



moment have the same overall geometry down locally to 
the level of [the coarser of] their cut-off lengthscales. For 
any given growing tip (labelled fc), the tip radius ak and 
flux density jk in the unprimed growth will be related to 



those in the primed growth by jfefl^ ^/^fc — j'k^"^ ^ 



where a is the local scaling exponent (as per Eq. (H)) of 
the harmonic measure between lengthscales and aj,. 
We take this exponent to have value a = atip on the 
grounds that this is locally a tip of the growth. Now let 
us focus on two different growing tips labelled 1 , 2 whose 
radii and flux densities are interrelated in the unprimed 
growth according to Eq. (||) by aij™ = 02 j^, and sim- 
ilarly in the primed growth by a[j'^ = fljj^™ ■ If we 
now insist that the advance velocities are in the same ra- 
tio (between tips 1 and 2) in both models, this requires 
{ill 32)^ — Ui/j'2)^ ' which forces the parameter relation 



l + m{l + a-d) l + m'{l + a-d) 



77 



rj' 



(4) 



For the two models to be equivalent in the relative ve- 
locities of all tips requires their parameters be related as 
above, where a = atip is the singularity exponent asso- 
ciated with growing tips. 

Although we have not strictly proved the equivalence 
of the models related above, we have shown that any 
such relationship must follow Eq. (^) and we will assume 
in the rest of this paper that this equivalence holds. All 
such models are then classifiable in terms of a convenient 
reference such as 770, the equivalent 77 when m — 0, cor- 
responding to the original Dielectric Breakdown Model. 
For example dendritic solidification with rj — 1 and 
777 = 1/2 corresponds to 770 = : it is thus not 

equivalent to DLA, but to another member of the DBM 
class. 
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Another puzzle resolved by our classification is a re- 
cent study showing conflicting scaling between DLA and 
different limits of a "Laplacian growth" model The 
latter model grows bumps of width proportional to flux 
density, so in the present terminology it corresponds to 
m = — 1 . The bumps are also grown with protrusion pro- 
portional to flux density. When the coverage of the grow- 
ing surface (per time step of growth) is low, then as the 
bumps are also distributed proportional to flux density, 
this limit corresponds to 77 = 3. By contrast high cover- 
age (with significant suppression of overlapping bumps) 
corresponds to 77 = 1. Using a = 0.7 (see below) these 
map through Eq. (|^) into 770 = 2.31 and 770 = 0.77 re- 
spectively, so the way their scaling was observed to 
bracket that of DLA is quite expected. 



III. THE ROLE OF NOISE 

DLA and DBM have been widely regarded as models 
in statistical physics, in that the local advance rate in 
Eq. has been implemented as the probability per unit 
time for the growth locally to make some unit of advance, 
entailing an inherent shot noise. Here we argue that dif- 
fusion controlled growth is a problem of turbulence type, 
with noise self-organising from minimal input. This was 
first suggested by Sander et al. Q but was only pursued 
in the case with surface tension cut-off, m = 1/2 in the 
present terminology, where it has been recognised more 
recently as chaotic viscous fingering p5| . 

The renormalisation of noise with lengthscale has hith- 
erto been discussed (at least for DLA) in the context of 
noise reduction |2^, |2^, focussing on the idea that 
as one goes up in lengthscale an equivalent coarse-scale 
model must have lower level of noise than crude shot 
noise. The limiting or "fixed-point" level of noise in DLA 
is small (at least according to references and ||l^) but 
certainly not zero, so it is natural to ask whether it can 
be approached from below as well as from above. The 
data in Fig. ^ show clearly that for DLA grown with 
very low noise by the methods of reference |13| , the rela- 
tive fiuctuations do indeed approach their limiting value 
from below as well as from above, and the same result 
was implicit in the earlier Renormalisation Group results 
of 11. 

The above discussion leads us to conjecture that for the 
full range of diffusion controlled growth under a contin- 
uum description of the interface, disorder in the initial 
conditions alone should suffice to feed instability, lead- 
ing to the same limiting levels of structural fluctuation 
on larger lengthscales as in the discrete models. The 
agreement we obtain below from simulations of the con- 
tinuum model without temporal noise provides direct ev- 
idence for this idea, which might also be argued obvious 
on the grounds that the MuUins-Sekerka instability |^] 
corresponds to a preponderance of positive Lyaponov ex- 
ponents in the dynamics. 
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FIG. 3; The size fluctuations ("output noise ) .Aout — 
{6N/N)'^, measured at fixed radius R for DLA clusters grown 
with off-lattice noise reduction at various ("input") noise 
levels Ain- For low Ain, Aout self organizes from below in the 
manner of a turbulent system. 



IV. CONTINUUM THEORETICAL 
DESCRIPTION AND SELF-ORGANISATION 

The new ideas above, that we can balance changing the 
cut-off exponent m by adjustment of 77, and that noise can 
be left to self organise, are the key to a new theoretical 
formulation of the problem, at least in two dimensions of 
space to which we now specialise. In two dimensions the 
Laplace equation in Eq. (^ can be solved in terms of a 
conformal transformation between the physical plane of 
z — x + iy and the plane of complex potential lu = (j) + i9, 
in which we take the growing interface to be mapped into 
the periodic interval 9 — [O,27r),0 = and the region 
outside of the growth mapped onto > 0. Then adapting 
reference Eq. (|^) leads to a closed equation for the 
dynamics of just the interface z{9) at (/) = 0, 



dz{9) 

at 



dz„ 
-'89'' 



89 
8z 



l+rj' 



(5) 



The linear operator V is most simply described 
in terms of Fourier transforms: 'P [J2k^~^''^ fk] — 

Efe P(fc)e-*'=«/fe = /o + 2Ef=i e-^'^Vfe, where we have 
introduced here an upper cut-off wavevector K. It is eas- 
ily shown that on scales of 9 greater than a smooth 
interface is linearly unstable with respect to corrugation 
for r] > (the MuUins-Sekerka instability ||l|), whereas 
for scales of 9 less than K^^ the equation drives smooth 
behaviour (corresponding locally to the case rj = — 1). 
This cutoff on a scale of 9, the cumulative integral of 
flux, corresponds in terms of tip radii and flux densities 
to aj « K^^, that is an m = 1 cutoff law. Thus the pa- 
rameter ?7 in Eq. (^ is more specifically 771 = ajyo, using 
Eq. (|) with d = 2. 

We will present a numerical study of Eq. (Bl) after a 



variable change in Sec. VII, where disorder was supplied 
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only through the initial conditions. The results clearly 
confirm that the equation self-organises critical scaling 
behaviour, without the supply of time dependent noise. 
The surprising form of the scaling is interpreted below. 



V. RENORMALISATION AND THE DBM 
SCALING LAW 

We now turn to a theoretical analysis of Eq. (||) , and 
for generality we will consider growth in a wedge of angle 
27rc (with periodic angular boundary conditions) so that 
c — > corresponds to growth along a channel whilst c = 1 
corresponds "radial growth" , that is growth out from a 
point in the plane. The primary theoretical requirement 
is that we must obtain results explicitly independent of 
the cut-off as A' ^ cx) : this is hard because we will see 
that the mean advance rate of the interface diverges as 
a power of K, and on fractal scaling grounds one would 
expect the same divergent factor to appear in the rate 
of change of other simple variables. One can of course 
take ratios of rates of change and look to order terms 
such that divergences cancel. To make this work we have 
been forced to introduce yet another change of variables, 

= ci?exp(ic6i- A(6i)) 



— cRexp ic9 — > A^e 



-ike 



fc>0 



which corresponds to Fourier decomposing the logarithm 
of the flux density. Here R is the effective radius of the 
growth and the non-analytic factor e^'^^ gives the mean 
winding of the conformal map through the wedge angle 
27rc, leaving X(9) as a simple Fourier series (except one- 
sided; see details in Appendix ^). The key to the success 
of the "logarithmic variables" A is that they decompose 
the flux density itself multiplicatively and, as we shall 
see, quite naturally capture its multifractal behaviour. 
In terms of these, time rescaled through dt = {cR)^ydt 
and y = (1 + r/)/2, the equation of motion (||) becomes 



^^-j:ik-j)Xk-,Pij)(ey('^-'^ 



+ 2 



(/fc-c) (e^(^+^))^ (7) 



where subscripts on bracketed expressions imply the tak- 
ing of a Fourier component, by analogy with A^. The ad- 
vance rate of the mean interface is correspondingly given 

by 



dR 
dt 



= cR 



(8) 



Details of the above analysis are given in Appendix 

At this point we can evaluate the multifractal spectrum 
in terms of these logarithmic variables. The multifractal 



spectrum of the harmonic measure follows from comput- 
ing the general moment Z{q,T) = J2 I'^^'l' li^-^r^, in 
the limit where all the intervals \6z\ and correspondingly 
66 approach zero; then the locus {q — l)D{q) = T{q) 
separates the limiting behaviour Z{q, t) oo from 
Z[q,T) — > 0. In our case it is convenient to fix 69 (ad- 
mitting wide variations in \6z\) and we must focus on 
the restricted range K^^ < 69 < 1. For 50 ~ 1 we 
have trivially Z ~ {cR)~'^ , whilst for 69 ~ the 
growth begins to look smooth so we can approximate 

Z{q,r) ElSer^l^r ^ K-^^^^^ Jd9\^\-\ 



Averaging gives 



dz I 



{cRY 



mg 



Z(g,T) ~if-«+^+i(ci?)-^(e 



„(A+A)t/2 



,(A+A)r/2 



yield- 



(9) 



The multifractal spectrum is then readily obtained from 
the separator behaviour discussed above, provided we can 
evaluate the average in Eq. (0). In Sec. VI we show how 



this can be done quite explicitly in the Gaussian Closure 
Approximation. 

We can now derive an elegant combination of Halsey's 
Electrostatic Scaling Law iQ combined with the tip scal- 
ing law of BaU and Witten ||, 0, that for DLA 
t(3) — Df = 1 -|- atip- This new derivation (unlike the 
earlier results) is not restricted to the m = case. The 
key idea is that we match the advance rate of the for- 
(6) ward tips of the growth (governed by atip) to that of the 
mean radius governed by a multifractal moment through 
Eq. (^). For consistency with the rest of this paper it is 
convenient to present the argument for the to = 1 repre- 
sentation, leading to tip velocity ^ = j'^^ ~ [K^^/aY^ 
where the tip radius a is set by the condition « aj ~ 
(a/i?)""p and hence ^ « i?-'?iif'?i(i/"tip-i). The over- 
all advance rate of the growth (in terms of its effective 
radius) is given from Eq. (|) by ^ « (ci?)"''! (e2'(^+^)). 
For the special value t/2 = y this can be substituted 
into Eq. (|), giving ^ w when Z is 

kept at its independent value (cR)^'^ . Comparing the 
two results leads to the DBM scaling law 



g(T = l + ?7i)=2 + ?7i/atip-2-h77o, (10) 
or its inverse form 



r(g = 2 + 770) = 1 + ?/i = 1 + atip??o ■ 



(11) 



The same result can be found, much more tortuously, 
from growth at general m. 

Now we turn back to the equation of motion of the 
logarithmic variables. Let us suppose some ignorance of 
the initial conditions and describe the system in terms 
of a joint probability distribution over the Xk, and let us 
denote averages over this [unknown] distribution by (..). 
We can in principle determine the distribution through 
its moments, whose evolution we now compute. For sim- 
plicity in this paper we assume translational invariance 
with respect to 9, so that only moments of zero total 
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wavevector need be considered, of which the lowest gives: 



-J2i''-j)Pij){Xk-,Xke^ 

j<k 

+ 2{k-c) (Xkef+^^)^+ic. conj.). 



(12) 

All of the higher moments lead to the same form of 

averages on the RHS, ^multinomial(A, A)e^(^"'"^^^, and 

all of these terms are conveniently expressed in terms of 
cumulants Q as detailed in Appendix ^ The key help- 
ful feature is that the expressions we require all naturally 



divide by one factor of ( e 



-y(^+^)'^ = ii(lni?), which is 

what we need in order to remove divergences by eliminat- 
ing increment of time dt in favour of id (In R) . The latter 
quantifies the update of large scale geometry, in terms of 
the advance in radius relative to the circumference of the 
wedge (or width of channel). The evolution of the second 
moments is then given in renormaliscd form by 



d{liiR) 



-^(A;-j)P(.7) (^Xk-jXke 
j<k L 

+ (c. conj.) 



-. „y(A+A) 



k 



(13) 



where ()^ denotes the Kubo cumulant p8|] . 

The above result (p^ ) is the key analytical step in 
this paper, because it removes divergent factors form the 
equations of motion. It is not dependent on the clo- 
sure approximation discussed below, and should support 
other approaches also. Moreover Eq. (p^), with the hier- 
archy of analogous equations for the evolution of higher 
moments, offers a new entry point towards exact results 
in the class of DLA and DBM models. 



where k* = c(l + l/?yi), and again the details are in 
Appendix ^ 

Equation ( p^ evolves to a unique steady state. The 
key to understanding this is to note that for k < k* 
the whole factor in large braces is negative definite, so 
S{k) = is the unique attractor. Then for k — fcmin, 
the first integer value above k* , S{k,^i^) = is unstable 
and the zero of the last factor leads to the global attractor 
having fcininS'(A:niin) = p-(fcniin — fc*)- The attractor valucs 
for higher k now follow by induction: denote the factor 
in large braces by B{k) and assume that the attractor 
has B(k) = and < kS{k) < p-, which are true for 
k = k^in. Then it follows that B{k + 1) = 2^ - kS{k) - 
{k + l)S{k + 1) and for kS{k) < the attractor must 
in turn have B{k + 1) = and hence a value of S{k + 1) 
conforming to < (fc -|- l)S{k -l- 1) < In the case 
kS{k) — ^ the only and stable solution is S{k -I- 1) = 
which leads to the same conclusions. Thus by induction 
the only stable attractor of the system has B{k) = for 
all k > k* and the corresponding steady state values are 



kS{k) = 0, 

m 

y2 

m 
y 



kS{k) 
kS{k) 



(^min ^ ) ; 

^{i + k*-k„ 



1 h ■ 

^1 "-mm 



3,. 



(15) 



Note that when fc* is integer, as particularly in the case 
c = corresponding to growth in a channel as discussed 
in fcmin = 1 -I- fc* and alternate valucs of kSik) are 
zero: in the case of channel growth this absence of even 
k is readily interpreted in terms of the dominance of one 
major finger and one major fjord. 

Within the Gaussian approximation and its predicted 
variances (15) we can now compute all [static] properties 
of diffusion controlled growth. From Eq. (^) we obtain 

Z{q,T) ~ (ci?)--exp(rV4Ef ^(fc)) ^ R-^k"^, 



using the values from Eq. (|15|), and hence Z(q,T 



8y^ 



R-''K 
haviour (now as K 



Thus the separator of limiting be- 
-> oo) is given by 



VI. GAUSSIAN CLOSURE APPROXIMATION 

To obtain simple tractable results we need to introduce 
some closure approximation(s) and we present here the 
simplest, neglecting all cumulants higher than the sec- 
ond, equivalent to assuming a joint Gaussian distribution 
(of zero mean) for A. This is entirely characterised by its 
second moments S{k) ~ (XkXk) which by Eq. (^ we find 
evolve according to 



_dS{k)_ 
'd(lni?) 



= 2y'Sik) I ^(fc - fc*) - fcS'(fc) 




q{r) 



1 



2(l + r?i)' 



(16) 



It is also easy to see that any closure scheme based on 
keeping cumulants of A up to some finite order leads to 
a corresponding degree polynomial truncation of ^(t). 

From the Legendre Transform of the inverse function 
T(g) (as detailed in Appendix |^) we obtain the corre- 
sponding spectrum of singularities, 



a I 









{v\ - 











(17) 



which in Fig. ^ is compared to measured data for DLA 
p9|, which later measurements l|30| reinforce. For the 
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FIG. 4: Multifractal spectra from the Gaussian theory (otip ~ 
2/3), compared to measured values for DLA Agreement 
is excellent for the active region t > 0, a < 1, and there are 
no adjustable parameters. 
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FIG. 5: Cumulative contribution to the mean growth velocity 
plotted against wavevector as fc''^ with logarithmic scales. The 
data are (bottom to top) for rji = 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 and 
all exhibit a common power law slope 1/a — 1 « 0.35±0.04 per 
the guidelines shown. The result were obtained by numerical 
integration of Eq. (flq) for c = (periodic strip geometry). 



region of active growth a < 1 {q > 0) the theory is quan- 
titatively accurate. At a = 1 it conforms to Makarov's 
theorernJM, and in contrast to the Screened Growth 
Model [|19|| it does this without adjustment. For a > 1 
the spectrum is only quahtatively the right shape, and for 
such screened regions our equations based on tip scaling 
may not hold. 

VII. SCALING PREDICTIONS FOR DBM 

To compare with the conventional DBM at to = and 
parameterised by 770 , we still need to compute theoret- 
ically the value of the tip scahng exponent which enters 
through 7ji = atipjyo- 

We can use the DBM scaling law ( p^ ) with the Gaus- 
sian Closure Approximation dig ) for (/(t) to fix the value 
of Qftip, and the resulting prediction is atip = 2/3 inde- 
pendent of rji. For DLA in two dimensions this value 
is respectably close to (but outside) measured values, 
atip = D — 1 = 0.71 ± 0.01 known from large direct sim- 
ulations of DLA ^ , but its suggested independence 
of 77 over a range of DBM is quite shocking. Numerical 
evidence, however, appears to lend support. 

We have investigated numerically what value of atip is 
seleced by the dynamics of Eq. (^ , with disorder supplied 
only through the initial condition. Changing variables to 

^ = (-zff (ci?)-ie— «)"''+"'^/' = ey\ we obtain 
dip dtj) d 

(18) 

where the rescaled time t is defined in terms of the evo- 
lution of cluster radius in Eq. (||). The tri-linear form of 
the RHS enables us to compute numerically the motion 
within a purely Fourier representation. 

In Sec. we have seen that the cutoff dependence 
of the tip velocity is u ~ i<f''i(i/"«p-i). This can be 



compared with the growth rate of the effective radius, 
^ (e!'(^+^)) = (V'V') = (V'V')o- So atip can be obtained 
from the if-dependence of (V'V')o (measured at fixed R). 
We can obtain it even from simulations with single K: 
the truncated sum Wcum(fc) = X]j<fe IV'il is expected to 
scale with k in the same way as the full sum {4'il')o does 
with K, because the Fourier components far below the 
cutoff should be insensitive to the value of K. Figure ^ 
shows the measured variation of Vcum{k) = X]j<fc I '^s 
k"^^ : this is expected to exhibit a power law with expo- 
nent (1/a — 1) and remarkably we obtain a « 0.74 ±0.02 
with no significant dependence on 771 in the range studied. 

It is a remarkable success for the Gaussian Theory to 
have predicted the completely unexpected insensitivity 
of atip to 7] . Whether this result can be truly an ex- 
act "superuniversality" is another matter, as certainly 
the Gaussian value for atip is only approximate and, as 
we discuss below, matters become more complicated for 
large 77. 



VIII. BREAKDOWN OF THE DBM MODEL AT 
LARGE 77. 

Sanchez and Sander |Q first noted that at high enough 
770 the DBM degenerates because all growth is domi- 
nated by the most active site, showing by direct simu- 
lations (at 777 = 0) that this happened around tjq w 4, 
a value reinforced by later discussions |3^ and new 
data [^. These discussions are particular to the to = 
representation and we believe can be associated with 
the degeneration of the moment governing the rate of 
gain of cluster mass: this scales with exponent r(77o) 
which degenerates to value 77oaniin when the moment 
becomes dominated by the (left) end point /(a) = 
of the multifractal spectrum. The fractal dimension 
given hy df — 1 + 77oatip — t(7;o) then degenerates to 
df = 1 -t- 7/0 (atip ~ ckmin)- If the least screened sites are 
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1 2 3 4 5 6 

FIG. 6: Singularity exponents reflecting the strength of 
screening as calculated from the Gaussian Closure Approx- 
imation with Otip set by the electrostatic scaling law. Note 
that for this theory Qtip > Qmin meaning that (within the 
theory!) the leading tips are not the most active sites. The 
screening transition arises when ascrocning governing the over- 
all screening hits ttmin, which it must subsequently follow. 
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FIG. 7: The fractal dimension for the m — Dielectric 
Breakdown Model as a function of 770 , using the Gaussian Clo- 
sure Approximation (curves) compared with published simu- 
lation data (points) |^^. The upper curve uses Qtip set by 
the Electrostatic ScaUng Law, which gives atip > Omin which 
is why df does not approach unity at the end point 770 ■ The 
lower curve shows how df does smoothly approach unity when 

we force Qtip = Qmin. 



the tips, atip = amin, then this also leads to d/ = 1 when 
V = Vo- 

The electrostatic scaling law leads to earlier transition 
in the behaviour, that is at lower 77, which is also more 
generic in that it does not depend on growth at some par- 
ticular value of m. This transition should also limit the 
applicability of calculating ^ exponents perturba- 
tively about rj^. The screening transition arises because 
the moment governing the mean screening of sites has ex- 
ponent T(2 + r]Q) which duly appears in the DBM version 
Eq. ( |ll| ) of the electrostatic scaling law, and this moment 
must hit the end of the f{a) spectrum before that cor- 
responding to t(77o) discussed above. Once we have hit 
this regime, at 770 > ?7o' ha,ve r(2-|-?7o) = (2 -I- 77o)ainin 
and the electrostatic scaling law degenerates to a form 
which can be rearranged to give 



1 



1 



1 

atip 



— • (19) 



This then leads us to choose between two scenarios: ei- 
ther (i) atip > Q!min in which case the behaviour remains 
non-trivial, or else (ii) amin = 1/2 which means the arms 
of the growth are essentially straight and we might sus- 
pect self-affine structure. 

The Gaussian Closure Approximation, with atip set by 
the electrostatic scaling law, has atip > cimin for almost 
all rj, leading to scenario (ii) above. Figure ^ shows the 
predicted variation of atip, amin and the value ascrocning 
corresponding to the exponent r(2 + rjo). The screening 
transition where ascrccning hits amin occurs (it can be 
checked exactly) at -qf — 2 corresponding to 770 = 3. 
Beyond this point atip is no longer quite constant and 
stays clear of amin, because of the changed functional 
form for r(2-|-77o), whilst of course ascreening follows a^^i^. 

The corresponding predicted behaviour of the fractal 
dimension calculated from df = 1 + 77oatip — T{r]o) is 



shown in Fig. ^ (upper curve) as a function of 770 up 
to 77q 5.4, where it has not fallen to unity because 
atip > ctmin is maintained. There is some change of 
functional form across 770 — Vo t>^^ i^ i^ scarcely no- 
ticeable graphically. Also shown for comparison (lower 
curve) is the behaviour when we force atip = Q^min in- 
stead of obeying the electrostatic scaling law: in this 
case the fractal dimension does smoothly approach unity 
as rjo rjQ, but unfortunately having sacrificed the elec- 
trostatic scaling law we cannot see anything relating to 
the screening transition. The predicted screening tran- 
sition at 77q = 3 seems to mark a break in the match 
to the simulation data of Hastings : below this con- 
forming to the electrostatic scaling law gives the better 
agreement, whereas beyond this better agreement comes 
from forcing atip — ctmin- The simplest interpretation is 
that the exact answer conforms to both conditions, and 
it is just their relative importance which changes around 
the screening transition. 



IX. THE PENETRATION DEPTH 

The multifractal spectrum suggests that the Gaussian 
approximation is good in the growth zone, so we have 
computed as a further test the relative penetration depth 
S, defined for DLA as the standard deviation ^ of radius 
of deposition divided by the effective radius R. For DBM 
more generally, we have for tractability used as measure 
the diffusion flux rather than the local growth rate. 

The key idea behind the calculation is that we cal- 
culate the relative distortion of the conformal map of 
the interface z{9) = i?e*°^ (l -|- X]fc>o wue^'^^) , away from 
the circular arc zq (9) = Re^"'^. Following Ref. p|, the 
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squared relative penetration depth is then given by 



i?2 



2tt 



de 



Re 



zi0)-zo{9) 



(20) 



,0 2tt \'^~ zoi9) 
which can be expressed in terms of the coefficients Wk as 



i?2 9 



Wk\ 



(21) 



fc>0 



From the definition 
identified as Wk = -z 



]) of X{9), the coefficients can be 
giving 



1 



fe>0 ^ ' 



de 

2^' 



ike 



d4) 

2^ 



(22) 



X e 



3 S p '^P 



This expression remains to be averaged over the distribu- 
tion of cluster geometries. Appendix ^ details the aver- 
aging of this under the Gaussian Closure Approximation, 
leading to the results shown in Fig. |[ 

For radial DLA the theory predicts S — (^/i?)r.m.s. = 
0.20 in rather modest agreement with 0.13 extrapolated 
from simulations 1 12 , and the prediction over a range of 



DBM parameter is shown in Fig.^. The predicted varia- 
tion of S for DLA grown in a wedge is also shown, and 
in the limit of zero wedge angle the value for the pene- 
tration depth relative to the width of the channel is the- 
oretically 0.13 compared with 0.14 measured [|2|. What 
is perhaps more interesting is the prediction of resonant 
features which arise when c, the wedge angle relative to 
27r, is an integer multiple of j^^, corresponding to inte- 
ger fc*, because the idea of resonant angles in DLA has 
been much discussed |l6[ |l^, this is the first direct 
and explicit prediction, which we look forward to seeing 
tested by simulations. 



X. GROWTH FLUCTUATIONS 



[I] 



2 
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FIG. 8: Upper panel shows the penetration depth relative 
to radius, for radial DBM growth as a function of r)i = a?7o- 
Lower panel shows the prediction for the DLA case, when 
growing in a wedge of angle 2-rc, plotted against c. The lower 
curve used r]\ = 2/3 and the upper curve used 771 = 0.71. The 
solid point corresponds to measurement for radial DLA, and 
the dashed line shows the limiting slope implied by measure- 
ments on DLA in a channel. Note the cusps predicted where 
c is an integer multiple of j+i^- 



of this under the Gaussian Closure Approximation is de- 
tailed in Appendix ^ 



Numerical confidence in the scaling properties of DLA 
was greatly bolstered by the idea of an intrinsic but low 
level of self-organised noise , so it is natural to ask if 
the present theory can address this. The simulation stud- 
ies of noise have rested on tracking the extremal radius, 
which is hard to extract from our analytic formulation, so 
we have had to compromise on something more accessible 
theoretically. The relative penetration depth has fluctu- 
ations which reflect the differing geometry of the growth, 
and for the case of growth in channel we have measured 
these fluctuations to be ((5S/S)siniuiation = 0.18. 

The corresponding theoretical calculation is a fairly 
straightforward generalisation of the penetration depth 
calculation itself: we simply calculate the average square 
of the expression in Eq. (p^, minus the square of 
Eq. (Dl) to obtain the variance of "E? . The evaluation 



The most useful comparison is for the channel, c — > 0, 
for which the penetration depth itself happened to be 
given rather accurately by the theory. In this case we 
obtain the relative variance of the square of the penetra- 
tion depth as {{£.IR)^) / {{i/Rff - 1 = 1.26, leading 
to (5S/S)thoory = 0.56 which is substantially higher than 
the simulation value. It is clear physically that the pene- 
tration depth comes predominantly from the lowest index 
modes of A, and this is apparent from our expressions 
above if we linearise Eq. (p^). Keeping only Aifor the 
channel would then make S,/ R be the magnitude of a sin- 
gle Gaussian distributed complex scalar, leading to very 
similar ^S/S. Thus it seems to be quite fundamentally 
the Gaussian form of our closure approximation which 
leads to overstatement of the penetration depth fluctua- 
tions. 
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XI. RELATION TO CONE ANGLE THEORY 



The angular resonances predicted in the penetration 
depth turn out to be in interesting correspondence with 
part of the earher Cone Angle Theory (CAT) of DLA 
||17| . In that theory a growing cluster was viewed as 
having an identifiable number of major arms n, and it 
was then further supposed that the growth should be 
marginally stable with respect to the loss of major arms 
through competition for growth. The strongest mode of 
such competition is where alternate fingers gain and lose, 
and the condition for this mode to marginally stable is 
in the present notation 
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n/c 



1 



1, 



(23) 



as calculated in ref for the case c = 1 and t^q = 1. 
In the CAT fractional n was presumed an acceptable ap- 
proximation and condition (|2^) was combined with ge- 
ometrical approximations to predict a, but here let us 
focus on the values of n implied. Using r\\ = arjo this 
gives 



n _ I 
2 " r?i 



(24) 



so our resonance condition corresponds directly to the 
case where the number of marginally stable major arms 
is an even number - which is of course required for the 
alternating mode stability calculation to be strictly ap- 
plicable. 

The CAT was closed in Rcf. by approximating 
the cluster as a solid polygon of n sides, for which a — 

leadmg to acAT = which is 

clearly quite different in principle from the Gaussian Clo- 
sure prediction of constant atip. However it is not easy to 
distinguish between them on the basis of previously pub- 
lished DBM data, as shown in Fig. H, and unlike CCA 
the CAT does not predict any other exponents. 



XII. DISCUSSION 

For DLA and its associated Dielectric Breakdown 
Models we have shown a theoretical framework which 
is complete in the sense that essentially all measurable 
quantities can be calculated. This extends to amplitude 
factors such as the relative penetration depth for which 
there is no theoretical precedent. For the full spectrum 
of exponents the practical advance over the Screened 
Growth Model is the elimination of fitting parameters, 
and it remains an open challenge to extend our theory 
to give quantitatively credible predictions for the large 
a part of the spectrum. For the exponent atip we have 
in the Gaussian approximation a striking new result that 
this is predicted constant over a wide range of rj, which 
begs direct confirmation by (expensive) particle-based 
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FIG. 9: The tip singularity exponent atip as a function of 770 
for the Dielectric Breakdown Model. Rising curve: Gaussian 
Closure Approximation; falling curve: Cone Angle Theory; 
points: simulation data for Qmin from reference [p3|. 



simulations. We look forward to addressing this in a 
following paper. 

For the DBM at high t] we find structure more rich than 
discussed hitherto, with a Screening Transition interven- 
ing before the upper critical value ryg is reached. Beyond 
the screening transition the scenario where a,nin = atip 
looks prospectively solvable (at least in terms of expo- 
nents) given the degenerate form of the Electrostatic 
Scaling Law which applies. The Gaussian Closure Ap- 
proximation leads to the richer but quite possibly mis- 
leading scenario amin < cttip, so sorting out the truth of 
this inequality would be very interesting. 

We have shown that the CCA naturally exhibits an- 
gular resonances which are in interesting correspondence 
with the ideas of the earlier Cone Angle Theory. No- 
tably these resonances now have clear predicted conse- 
quences such as we demonstrated for the relative pene- 
tration depth, and they can be explored either by grow- 
ing in a wedge of variable angle or by varying rj - so once 
again behaviour vs 77 is a key probe of our understanding 
of the problem. For DLA in particular the best theoret- 
ical value of atip remains 1/^/2 « 0.71 from the CAT 
Ip^l , but the greater test now lies in the 77 dependence on 
which CAT and GCA differ qualitatively. 

Within DLA and DBM we look forward to calculat- 
ing more properties such as the response to anisotropy, 
which is fairly readily incorporated into our equations of 
motion. The hard part is that in breaking angular sym- 
metry we can no longer exclude non-zero first cumulants 
(Afe), and a full matrix of second cumulants, but the cal- 
culation is in principle straightforward. A conceptually 
more challenging avenue is to improve on the Gaussian 
approximation itself which we have used to obtain ex- 
plicit theoretical results. Truncating at a cumulant of 
higher order than the second is hard, and more seriously 
it does not correspond to a positive (semi-)definite prob- 
ability distribution. An alternative route of improvement 
which we are exploring is closure at the level of the full 
multifractal spectrum. 
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Whilst our main use of the equivalences within the 
class of (?7, to) models has been to facilitate calculation, 
through mapping onto m — I, the particular associated 
claim that surface tension control is included through 
TO = 1 may prove controversial. This would imply that 
the scaling properties of the chaotic viscous fingering 
regime can be predicted from suitable DBM simulations. 
The DBM simulations required are relatively accessible 
and the greater difficulty in pursuing this agenda lies in 
obtaining suitably calibrated experimental data or accu- 
rate direct simulations of fingering out to high degrees of 
ramification. 

There are possibilities for wider application of ideas in 
this paper, where we have formulated DLA and DBM as 
a turbulent dynamics governed by a complex scalar field 
in 1+1 dimensions. Decomposing this field multiplica- 
tively (through Fourier representation of its logarithm) 
was the crucial step to obtain renormalisable equations 
and theoretical access to the multifractal behaviour, even 
though other representations offered equations of motion 
( p^ ) with weaker non-linearity. It is natural to specu- 
late whether the same strategy might apply to turbulent 
problems more widely, where the key issue appears to be 
identifying suitable fields to decompose multiplicatively 
which are of local physical significance, and subject to 
closed equations of motion. 



2,3,..., the prefactor in front of e~'^" has to vanish. The 
Gaussian Closure solutions (p^ ) automatically satisfy this 
constraint, because they have Afc = for fc < c. 

In the region far from the growth the leading term 
dominates, giving z{uj) = Re'^'^ which shows the sig- 
nificance of R: it is the apparent radius of the growth 
(corresponding to (/) = 0) as seen from far away. 

To obtain the transformed equation of motion (|^) , first 
take the logarithm of Eq. (||), giving —X{0)+hiR = 
In (ll) -I- In (— i/ce~*'^^). Then differentiating both sides 
with respect to time (at constant 9) gives 



dX 
~dt 



OR 
Rdt 



dz\ ^ d dz 
m) de'dt 



{cR)-^yeyi^+^) 

(ci?)-2i'e^(^+^) 



(A3) 



where we have used \d9/dz\ = {cR)~^e^'^ and subse- 
quently the powers of R can be taken outside V. It is 
then trivial to take Fourier components of both sides to 
obtain Eqs. (|^) and (^, the latter coming from the zeroth 
component which we chose to be absent from A. 
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APPENDIX A: LOGARITHMIC VARIABLES 

In this Appendix we present some details about the 
logarithmic variables in Eq. (^ and their equation of mo- 
tion (|) and (|). 

We begin by explaining the choice of the analytic form 
in Eq. (ph. In terms of w = (fi + id, in which complex 
plane the region exterior to the growth is mapped to a 
[half] strip, Eq. (g) analytically continues to 



— = cRe''^e-^''>o: Mc^*" 
duj 



Afce 



(Al) 



Expanding the second exponential factor above to all or- 
ders gives 



4^ = ci?e^"(l 



Ai e 



(A2 - A^/2)e-' 
(A3 - ...)e-3" 



(A2) 



and integrating with respect to tu gives a conformal map 
from the right hand part of the strip to the exterior region 
of the wedge as required. When c is positive integer, 
periodic boundary conditions put a restriction on the A's: 
Al = when c = 1, and for the less interesting cases c = 



The Kubo Cumulants for independent variables Xi 
are given in terms of their moments by In (e^» = 
^e^i — l)^ , where Pi are arbitrary (scalar) pa- 
rameters |2^. From this it is well known that 
by differentiation with respect to parameters one ob- 
tains (Xie'^) = {Xie"^) (e^) and (XiXae^) = 

((M2e^)^-f(Xie^)^(/2e^)J(e>^), where W ee 
PiXi, and similar results for higher moments such as 

To obtain the renormalised equation of motion 
( [Tsl ) from the unrenormalised equation ( p^ we 

need first to apply the above to (^■^ke-^!'''^^^^'^ ■ 

Writing this as / ^e**^^ <^Afcey(^('')+^(''))^ allows 
us to apply the cumulant identities directly, giv- 



ing 



-j 



In the trans- 



hence J2j \^kel 

lationally invariant case considered here only the 
J = term survives in the latter summation, lead- 

ing to (A.ef+^)) . (A.ef+^))^(er^"^) as 

required. The calculation of (|Afc_j AfeeJ^'^^^''^ = 

A.-,A.ef +^))^ + {x,-,eff'^)^ (A.ef +^)) J x 
vi^+^)\ pi-oceeds along precisely analogous fines. 
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APPENDIX C: LEGEND RE TRANSFORM OF 
AN INVERSE FUNCTION 

In Eq. (|l^) we are confronted with a slightly unusual 
situation: we wish to find the Legendre Transform /(a) 
of the function T{q), that is f{a) — qa T{q) where 
a = dr/dq, given a simple form for the inverse function 
^(t). Let g{x) be the Legendre Transform of q{T), that 
is g{x) — XT — q{T) where x = dq/dr. Thus x = l/a 
and we have g{l/a) — r/a — q — —f{a)/a, so the two 
Legendre Transforms are very simply related. 

The example needed from Eq. ( ]l6| ) has the form g(r) = 
1 + T + br^ leading to x = 1 + 2&t and hence g{x) — 



-1 + bi 



2b 



. Then we have /(a) = ~ag{l/a) leading 



directly to Eq. (17). 



APPENDIX D: CALCULATION OF THE 
PENETRATION DEPTH 



We start from the expression given in Eq. ( |22|) for the 
square of the relative penetration depth without any av- 
eraging over clusters. As the expression is simple expo- 
nential in the it is straightforward to average over a 
Gaussian distribution leading to 



^2 



' k>Q ^ ' 



2- da 



Note that the term —1 in the integrand makes no difference for fc 7^ 0, but we can include fc = in the outer 
summation. This enables us to rearrange the combination of summation over fc and integration, giving 

S2 ^^J^°°dxxe^^ (^e^.sip)e-''^ _ ^ ^2) 
as is readily verified upon expanding the exponential e^p ■^(p)^ f^Q gXl orders. The summation in this exponential can 



be evaluated in closed form using the variances from Eq. ( |15D and the standard forms u + u /3 + u'^/5 + = In y 

and u'^/2 + u'^/4: + u^/6 + .... = In y/j^r . This leaves one numerical quadrature to obtain the results shown in Fig. ||. 

To compute the fluctuations in relative penetration depth we again start from Eq. ( p^ ) and now average its square 
to give 

(D3) 



\^«/^^ / 4^2^ (c^kf J Jq 27r27r ^ ic-k'f J J. 27r 27r 

X ( exp ( - ^ \pe-'P^ - ^ Ape'P* - ^ A^e" 



p 




where averaging the last factor gives exp {y,^ S'(p)[e^P('^~'') + ^ip(4>-e+0-e') ^ ^zp(4>'-e'+0'-0) ^ g«p(</.'-e')]^ ^ q 

tegration with respect to an absolute angle is now redundant, and the integrations with respect to (f) — 9 and 0' 
can be rearranged using the same trick as in calculating ((^/i?)^) previously, leading after some cancellations to 



nc m- 




{D4) 



dx'x'e''''' I — LY.^S{p)[e-''-+e-'"^--^^>+e-''^-'+^^>+e-''-'] _ gEp S(p)[e-P-+e-''==' ] 

2n V 



Above we used the simplification Jq'^ ^e^''^^^'''^ ^ = 1, which follows upon expanding the exponential 

where only the leading term survives. This can be further applied leading to the more compact form (5(2^)^ = 



S(p)e- 



1 



Particularly in order to address the limit of channel 



growth, c — !■ 0, it is convenient to bypass normalisation conventions by looking at the relative fluctuations in ^ 



"2 
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which are now given by 



■=•2 



S(p)e-P== 



(D5) 



From this we obtained the resuh cited in Sec. |x| by numeri cal quadrature. 
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